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NONSmTIOHARY GAS FLOW IN THIN PUTS 
OF VARIABLE CROSS SECTION* 

By G. Guderley 


ABSTRACT 


Characteristic methods for nonstationary flows have "been 
published only for the special case of the isentropic flow up 
until the present, although they are applicable in various 
places to more difficult questions, too . The present report 
derives the characteristic method for the flows which depend 
only on the position coordinates and the time. At the same time 
the treatment of compression shocks is shown. To simplify the 
application numerous examples are worked out. 


*'Niohtstationare Gasstr'dmungen in dunnen Rohren veranderlichen 
Querschnitts • " Zentrale fur wlssenschaf tliches Berichtswesen der 
Luftfahrtforschung des Generalluftzeugmeisters (7MB) Berlin -Adlershof , 
Fore chungsber ioht Nr. 1744, Braunschweig, Oct. 22, 1942 
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1. INTRODUCTION 


In papers by F. S chult z-Grunow and R. Sauer methods have 
been developed recently for completely solving the problem of 
nonstationary isentropic gas flows in a pipe of constant cross 
section. An expanded view of the problem is the basis for the 
present report. Flows are considered, which likewise depend only 
on the position coordinate; however, the crosB section of the tube 
need no longer be constant and the entropy may vary from particle to 
particle. The method of solution applied here has been discovered 
almost simultaneously in several places, by Adam Schmidt, W, During, 
and F. Pfeiffer, among others. 

The application of the characteristic method is possible 
without a previous substantial knowledge of mathematics. Corre- 
spondingly, if a derivation was desired too, one could be had 
which did not make any special mathematical demands on the reader. 

As a model, the Busemann derivation of the characteristic method 

for two-dimensional stationary gas flows might possibly do . It is 
actually possible to apply this derivation immediately to the 


Schultz -Grunow, F. : Nichtstationare eindimensionale Gas- 

bewegung. Forschung auf dem Gebiet des Ingenieurwesen^ Bd. 13' 

(1942) pp. 125 to 134. Sauer, R.: Charakteristlkenverfahren fur die 

e indimens ionale lnstationare Gasstromung, Ingenieur-Archiv, 

XIII Vol. (1942) pp. 79 to 89. Vorbereitende Untersuchungen sowie 
Anwendungen finden sich in den Arbeiten von H. Pfriem. Zur Theorie 
ebener Druckwellen mit steiler Front Akustische Zeitschrift 
Jahrg. 6 (l94l) part 4. — Die ebene ungedampfte Druckwelle grosser 
Schwingungsweite, Forschung Vol. 12 (194) pp. 51 to 64 — 
Reflexionsgesetze fur ebene Druckwellen grosser Schwingungsweite, 
Forschung Vol. 12 (l94l) pp. 244 to 256 — Zur gegenseitigen 
Uberlagerung ungedampfter ebener Gasswellen grosser Schwingungsweite, 
Akustische Zeitschrift Jahrg. 7 (1942) part 2 — Zur Frage der 
oberen Grenze von Geschoesgeschwindigkeiten Zeitschrift f. techn. 
Physik 22 (l94l) pp. 255 to 260; 'Eine weitere Anwendung findet 
sich bei G. Damkohler und A. Schmidt, Gasdynamische Beitrage zur 
Auswertung von Flammenversuchen in Rohrs tr.ecken. Zeitschrift 
fur Elektrotechnik Vol. 47 (l94l) pp. 547 to 587. 

^Busemann, A.: Beitrag Gasdynamik in Handbuch der Experimental- 

physik (Wien-Harms) Bd, 4, Tail 1, _p. 421 and adjoining pages. 
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isen tropic nonstationary flow in a pipe of constant cross section 
and from this by means of some supplementary physical concepts 
succeed in getting a treatment of flows in a tube of variable cross 
section; this is the course which had been taken, originally. Jn 
comparison to the mathematical theory of characteristics , however, 
these considerations operate with a lack of clarity sufficient so 
that the mathematical theory - for the engineer too - can be 
represented as the best approach to the characteristics method. 

Considerations necessary for the present problem are now 
brought forward frcm the characteristics theory 3 . As a result, 
equations for the directions of the characteristics as well as 
conditions which must be satisfied along the characteristics are 
obtained. Proceeding frcm these relations, the next sections 
develop the actual method of computation. Next, the characteristics 
method for the case which is familiar by now, that of isentropic 
flows in a pipe of constant cross section, is deduced again and 
the transformations appearing there are used to simplify the 
computation in complicated cases, too. Since this is not always 
possible, the most general form of the characteristics method is 
shown in a later section. After this, the formulas obtained for 
the special case of an ideal gas with constant specific heat' are 
simplified and the consideration of boundary conditions explained. 
The remaining sections deal with calculation of compression shocks; 
the known relations which connect the phase before and behind a 
compression shock with one another are set forth in a convenient' 
form for the present problems and with that the calculation of a 
compression shock in a flow is carried out . 

The theory is illustrated with suitable examples treated in 
detail . In that regard. It seemed advantageous to avoid definite 
problems of technical interest, in doing so gaining the possibility 
of working out examples under very general assumptions without 
excessive effort. It Is hoped that, nevertheless, the application 
of the method to physical problems offers no additional difficulties 
worth mentioning inasmuch as the earlier publications contain such 
applications. The author expresses his thanks to Dr. Fa ne Lehmann 
for working out the examples. 


3 11 f " t " 

Compare Courant— Hilbert, "Methoden der mathematischen Physik II", 
p 291, Guderley follows the representation given by H, Seifert at the 
same institute for Gas Dynamios in lectures. 
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2 . BASIC EQUATIONS 


Consider nonstationary, perfect gas flows in a pipe with a 
cross-section that varies in space and time 4 in the neighborhood 
of the flow tube; that is, it is assumed that the velocity and 
the phase over a cross section of the pipe may be considered as 
sufficiently constant. In general, this assumption is Justifiable 
only if the thickness of the tube, relative to its length, changes 
slowly enough . Only for flows which have as surfaces of constant 
phase parallel planes, coaxial cylinders or concentric spheres 
need this limitation be ignored. These flows with plane, cylindrical, 
or spherical wave propagation are included as special cases in the 
present problem. 

To stress the relationship to stationary two- dimensional flows, 
let the axis of the pipe be vertical, the position coordinate be y, 
the time be t and plotted horizontally. In this yt- diagram the 
flows are investigated. (Compare fig. 1 .) 

Let 

p pressure 

s entropy per unit mass 

p density 

v velocity 

F a F(yt) the cross Beotian of the pipe, let F be given 

In a region free of compression shocks, the flow is described by 
the dependency of the density on pressure and entropy, the Newtonian 
principle, the equation of continuity and the statement that the 
entropy of a particle is preserved, as follows: 


P = p(s,p) 

(1) 

l&U T gz + |j» 0 

p oy dy dt 

(2) 

v|£ + p|z + |£+ P vto+ p 0 

dy dy dt dy ^ St 

(3) 

„ Sb . Sb a 
5y St 0 

(4) 


Problems with 'time variations in the cross-sectional area 
are rare* they were included, since they can be handled without 
additional difficulty. 
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In this the derivatives of p should he replaced hy the derivatives 
of p and s, for this purpose 


dp 1 
3 * ‘ a 2 


(5) 


is introduced. Therefore, instead of equation ( 3 ) 


X.^ + T | £ ^. + | S+ l^ + ^|a + aM 

a 2 dy ds dy dy a 2 dt ds dt dy 


dlnF , din I ~ / Oo \ 

-5=- + P “gjr * 0 (3a) 


is obtained. 


3 . FROM THE THEORY Of CHARACTERISTICS 


In regard to the system of equations (2), (3a), and (^), the 

familiar question is raised from the theory of characteristics . In 
a region of the yt~plane let the solution of this system of 
equations and its derivations be finite throughout . On a curve C 
placed in this range let the values p, v, and s which correspond 
to this solution and, therefore, the appropriate derivatives taken 
in the direction of C be known. The question is asked whether the 
derivatives in other directions may be computed with the aid of the 
system of differential equations, and under what conditions. To 
answer this, a curvilinear coordinate system is introduced 
in which a curve £ = constant coincides with C (fig. 1). All .the 
derivatives with respect to ^ along this curve £> -=. constant are 
given, -the derivatives with respect to are sought. This trans- 
formation Is carried out and terms are arranged so that the unkn own 
derivatives with respect to are on the left and only known 
quantities are cai the right. That Is 

S - 

^ = T)(y,t) 

. . £>_.=, 3§ + A A). 

dy dg dy dTi dy 

+ A 'da. 

dt d& dt dji dt 
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With this form (2), (3(a)), and (!+) are obtained 


+ + §£\ = _^!£ i _ $l( r ?£l + £ll\ 

dg p dy 5g \ dy 5t / 5 tj p Sy dp V dy dty 


dp 1 

dla 2 



as 

Sy 


4- 




£ 

Sy 


as $e/t $1 + &l\ 

dg ds\ dy dty 



dp \ _Sv dp _ bal dp , dp\ dp _ 
dt/ dp P 5y dpy dy btj ds 


pv 


bln F dZnff > 

dy P St [ 


2§f T 11+ 21) . J» ( T 2a + 22 ') 

d£ dy Sty S’! Sy Sty 


A linear system of equations is obtained for the unknowns ~ , 

5s 

and the unknowns, themselves, are obtained by Cramer's rule as 


the quotient of two determinants . The determinant in the denominator 
is the same for all unknowns . It always gives a single -valued 
so3.utian for the system of equations, if the determinant in the 
denominator is different from 0. In the other case with a determinant 
in the denominator that vanishes, it is a necessary condition for Hie 
existence of solutions that remain finite that the determinants in 
the numerator also vanish. In this case, however, the solution of 
the system of equations is only defined over any portion of the 
solution of the homogenous system. In application to the By stem of 
equations (6) signifies the following: The determinant in the 

denominator is formed from the coefficients of the unknowns. 

Considering a fixed point on C, at which p, v. and s are 
known by assumption, the coefficients depend on £b dg that 

dy dt 

is the direction of C . If C Is so directed that the determinant 


in the denominator does not vanish anywhere, the 


etc . are 

d£ 


computable as single valued. 


Of greater interest for our considerations (a the other case, 
namely, that the determinant in the denominator Is zero at every 
point of C . Such a curve is termed a characteristic . Because of 
the assumption of finite derivatives the determinants in the 



NACA TM So. 1196 


7 


numerator also vanish. Relations are obtained thereby, in which the 
right side of (6) and, therefore, the derivatives of p, v, and s 
along C appear as essential ingredients. These relations represent 
the starting point of Idle graphical numerical method of solution. 

Since the solutions of a linear system of equations are no 
longer single valued for vanishing numerator and denominator 
determinants, the pursuit of a given solution of a characteristic 
is possible in various ways. These different possibilities actually 
appear on changing the initial and boundary conditions. 


k. THE DIRECTIONS OF THE CHARACTERISTICS 


To find the directions for which the curve £ » constant is a 
characteristic, the determinant in the denominator must be set equal 
to zero in the solutions of the system of equations (6) . 



Prom this are obtained -the conditions 

T ! + !t " 0 (7 > 


or 


L + »\ * 31 . 0 

V 75 S ° 


(8a) 


I 
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or 


(v - a) 


ag + as 


o 


The slope of any curve g = constant is given hy 


(8b) 


dy _ 8t 
. dt 

dy ' 


From. (7) and (8), together -with this, the slopes of the charac- 
teristics are 

& = v (9) 

df 

or 

dy. = v + a (10a) 

dt 

or 


& = v - a (15b) 


The characteristics defined by (9) are path- time curves for the 
individual gas particle) they might be termed life lines of tho 
particles. According to (10), velocities are determined from the 
slope of the other characteristics, which differ from the velocity 
of the particles by ±a- For stationary flows the Mach waves 
correspond to these last characteristics) this designation will be 
adopted. Therefore, let Mach waves of the first family be those 
which spread out with the velocity v + a and Mach waves of the 
second family be associated with the velocity' v - a. 


5. THE COTSISTMCY CCHDITICNS CN THE CHAMCTERISTICS 


As shown in section 3 , along the characteristics, certain 
conditions must be complied with by the derivatives which result 
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frcm the vanishing of the determinant in the numerator. These 
conditions are called consistency conditions, for p, v, and s 
are subject to them, if the derivatives with respect to £ are 
to remain finite. If the right side of (6) is designated B 1# Eg, 

and Bo in sequence, then the following is obtained for the 
determinant in the numerator of the quotient for 2E: 


\ 

T l§ + P 

Sy ot 

0 



*2 

dy 

|e. (v 

os dy dy 

= 0 

(11) 

*3 

0 

f T |i + |iV 

\ oy 5t J j 




This determinant must vanish to give the directions of the charac- 
teristics. Substituting (7) gives 


E. 


E r 


E. 


0 

a| 

5y 

0 


0 

0 


= 0 


According to this the determinant (ll) vanishes by itself. With 
(8a) , that Is, for a Mach wave 1 



is obtained, or 
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In this, S. is certainly different from zero, as long as v and a 

dy 

are finite and grad £ y» 0 . As the condition for Che Mach wave 1 is 
obtained 

-pE 1 - aBg + a H3- « 0 (12a) 

The consistency conditions for the Mach waves 2 is, if a is 
replaced by -a 


do 


-pE-l + aE 2 - a U E 3 = 0 


(12b) 


A condition for the life line is obtained if the vanishing of the 
determinant in the numerator in the .'Quotient for ^s to be got 

from (6) requires 


1 dt 

(a£ „ aS\ 

P dy 

n N ?r + Sty 

v ^ . e£\ 


v dJ 5t) 

0 

0 ¥ 


B, 


E, 


E. 


= 0 


For the Mach wave this equation is satisfied by itself, the condition 
for the life line is 


b 3 = 0 


(13) 


The determinant in the numerator of could be investigated, too; 

d& 

however, this would not give any new consistency condi ticns. 
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t 


The values of R^, 
Frau (12a), 


R 


2 , and R^ are still to he put in. 


1 dp 

a 5? 


|a . a) 


+ la 

bt 


hr 


^3. (v + a) + j s -apv S^B? 


dlnF _ „ dlnF 

— a -sr 


is obtained. The direction 
(10a) j on. that account 


along a Mach wave 
dt 


1 is given by 


da _ Sa djr 
dt dy dt 




is valid for it. With that the consistency condition for the Mach 
wave 1 can be written in the form. 


1 dp + dv 
ap dt dt 




(14a) 


The consistency condition for a Mach wave 2 is obtained, by 
substituting -a for a 


2_ §£ 

ap dt 


dv _ 
dt 




From (13) for a life line is obtained 


(ito>) 



This may be integrated immediately 


s = constant ( 15 ) 

Naturally this constant will differ from, particle to particle, in 
general . 



12 


WACA- m Ho . 1196 


6 . FLOWS WITH CCN STMT 0TTROPY ML IK IFORM PIPE CROSS SECTION 


Equations (14), (15), (9), and (10) just obtained are certainly 
useful, fundamentally, as a starting point of a characteristics 
method - in fact, there are "examples, where it is necessary to revert 
to them (compare section 11); in most cases, however, there are 
still other transf ormations suitable . The direction in which to 
proceed for these iB obtained if an attempt is made to derive the 
characteristic method for isen tropic flows in a pipe of uniform 
cross Bection from equations (14) and (15) possibly in the form 
applied by Schultz -Grunow . _ To emphasize the fundamental ideas, no 
assumptions of any kind are made therein of the characteristics of 
the flow medium. 

On account of the hypothesis of constant entropy, equation (15) 
satisfies itself . In equations (l4) the right sides are omitted 
since it concerns a tube of uniform crosB section. Further, on 
account of the hypothesis of constant entropy the state of the gas 
is still dependent as only one variable, perhaps the pressure, 
or the temperature } the quantities appearing on the left side 
and a are accordingly functions of this variable. It is possible, 

therefore, to consider the expression dp. as a differential. Let 

pa 


T temperature 
i heat content (enthalpy) 
s entropy 

By the second main theorem 

T ds = di - i. dp 
P 

From this, on account of the hypothesis of constant entropy, 



With that, it follows that 


dp 

pa 



dT 
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w(t) 


JL / 1 ( ^£\ cLT 

a o J? aV '^s 


(IT) 


is introduced in which. a Q is the sonic Telocity of a comparison 

phase ■which was added to make W dimensionless . The phase of the 
gas may he characterized hy W from W = W(T). It follows' that 

T = T(W) (18a) 

Further, it is valid that 

p = p(T) = p(W) (l8h) 

a = a(T) = a(W) etc. (l8c) 

With the use of W equations (li)-) appear in the form 


* 

a Q dW + 

dv = 0 

for 


a Q dW - 

S 

II 

o 

for 


Bringing in 

X 

= W + 




= W 


V 


(19a) 


(19b) 


these last relations change to a form which may he integrated. This 
gives 


X = constant 


for Mach wave 1 


(20a) 


\x ~ constant 


for Mach ware 2 


(20b) 


If the magnitudes of X and n are known for a point of the yt-diagram, 
the velocity is thereby completely defined as well as the thermodynamic 
phase. It is, to be exact, 


W = LJUi- 
2 


(21a) 


X = l^ZJL (21b) 

a o * » 
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and on account of equations (l 8 ) 

P = p(X + u) (22a) 

a = a(X + |i). ( 22 b) 

The next sections explain this transformation and the application 
of equations (20) to an example of an. ideal gaB whose specific heat 
is a function of temperature . 


7 . THERMODYNAMIC RELATIONS FOR AN IDEAL GAS WHOSE SPECIFIC 
HEAT IS A FUNCTION OF TMPERATURE; COMPUTATION OF W 


Let 

c^ specific heat at constant pressure 

c T specific heat at constant volume 

R gas constant. 

For an ideal gas 


£ = RT (23) 

P 

According to the second main theorem, if p and T are considered 
as independent variables 


ds = 1 dT + i. c|i_ <jp m 1 ip 
T ST T op pT 


(24) 


Since ds is a perfect differential. 


8 (l di\ 

5? \t St/ 


S3 


5 ( 1 8 i 
§T \T 5? 



Accordingly, substituting p from ( 23 ), the following known fact' 
is obtained 
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1 


that is 


now 


therefore 


From (2h- ) as a result 


i . i(T) 



c (T) 
P 


ds = JL cLT r R 4 d. 


T P 


• (25) 

(26) 


(24a) 


and from this by integration 




22. 
. R 


-In JL 
T Po 


the index o characterizes a comparison phase • 


Introducing 


.T 
■ > 



2p ar 
r t • 


(27a) 


n 



(27b) 


gives 


-E- = Jit (28) 

Po 


Considering c.p(T) as knovm, p by (28) and p by (23), 

are given as functions of T and s ; the thermodynamic properties 
of the medium can be calculated in principal, therefore. The 
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quantity a is also defined by p and p . The computation of a 
is by all means simpler if carried Cut in Hie following way. 
According to (5) 



for which the entropy is to be kept constant. For constant entropy 
from (24a) 


dT _ _R 
T " c p p 


by differentiation from (23) 


dp _ dp dT 
P " P = T 

From the last two relations together with the familiar relation 


is obtained 


Cy “ F 


a(T) =/S 2. ET (29) 

' c v 

The relations discovered up until now describe the properties 
of the gas and must always be known 5 it makes no difference which 
variation of the characteristic method is chosen for the calculation 
'of the flow. In contrast, the introduction of the functions W, X, 
and p serve only as preparation for carrying out of Hie charac- 
teristics method in the form presented in the preceding section. 
Next, to compute W. From (25) it follows 



f|i\ - 

\3T L 


with (26) 
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Setting the last equation as 


\f — 1 + u 


For the present case, s = constant = s Q (28) becomes 

'.T . ' 

dT 

T R T 

JEL (T) = e (28a) 

P ° 

Wow the following can be formed , • 

T = T(X + p) 

a * a(X + p) 

P = P(X + u) 

These calculations were carried out numerically , for carbon 
dioxide. The relation between the specific heat and temperature was 

taken from Eiitte with the aid of these values (i - i Q ) /a D ^ , 
a/a D , P, and W can be computed from equations (26), (25?), (27a), 
and (30) as functions of the temperature. (See figs. 2(a) and 2(b).) 
Figure 3 shows a/a Q , P, and T plotted as functions of 
X + p = 2W. 

8 . TEE! CONSTRUCTION OF TEE FLOW FIELD 



well as (29) into (17) gives 


-i/V 

a o Jt„ Y 


c P c v 

RT 


dT 


(30) 


The following problem should be dealt with: Along a curve K 

of the yt-diagrara, which has at the most one point in common with 
each characteristic, let p/p Q and v/a 0 be given (fig. 4) . The 
flow should be constructed for the following times as far as it is 
defined by the portion of K given. Therefore, it is concerned 
here with the computation of the part of the flow defined by the 
initial conditions which by the same arguments appear everywhere .in 
the int erior, too. Before the construction of the flow can be started 

8 

fiutte, 27th edition, Vol; 1, p. 48, table 5, Berlin 1941, 
Wilhelm Ernst und Sohn, publishers . 
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the initial values p/p Q and. v/a Q must be expressed in terms of 

the variables X and (i • Since the entropy was assumed constant, 
p/p o and P coincide . With the aid of the relation presented in 

figure 3, between P and X + p and eguation (21b), X and p 

may be ascertained without difficulties . Figure 4 shows on the right, 
the yt-diagram, on the left, the diagram of the assumed values p/p 0 

and v/a Q as well as those of the computed quantities X and 4 as 
functions of y. 

Proceeding from the individual points of K if the network of 
Mach waves had been brought in, the phase at each lattice point 
would be determined thereby; according to (20) X is constant along 
Mach wave 1, p along Mach wave 2, and on that account, equal to 
the values at those points of K from which the Mach waves spread 
out. By (21b) and (22) the phase is given by X and p. To be 
able to draw the network of Mach waves, only their directions are 
still needed. These are given at the lattice points by (lO)j 
a/a 0 is a function of X + |i in figure 3, v/a D is computed 



The direction for the portion of a Mach wave between two lattice 
points is approximated as the average value of the corresponding 
directions at the lattice points . 

The construction becomes especially simple if the Mach waves 
are drawn for equidistant values of X and p. The directions 
of the Mach waves appearing can be computed beforehand and possibly 
prepared in the form of table I. The interval between adjacent 
valueB of X or p was selected as 0 .1, the size of the interval 
depends cn the accuracy desired. In the table the upper column 
headings and signs refer to Mach wave 1, the lower to Mach wave 2 . 

The numbers entered in the table represent the average values for 
(v + a) /a and (v - a)/a Q . For Mach wave 1 for which X = 0.3 
and which leads from a point with p = 0 .2 to a point with p s* 0.1, 
in the column with the heading 1 = 0-3 the value is to be taken 
frem the row p = 1.5, that is, (v + a)/a Q = 1.103* 

In the flow diagram, the values of X valid there are entered 
to the left of the lattice point and the values of p to the 
right. To determine, for example, the position of C from the 
■points A and B since the phase of C is given beforehand by 
X = 1.1 and p = 0.5 the average directions of the Mach waves 
(v + a)/a D = 1.422, (v - a) /a Q = -0.778 can be taken from table I 

and drawn in the yt-diagram. The auxiliary diagram on the left in 
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figure U- can "be used for this. There the direction of a Mach wave 
for vhich (v + a) /a 0 =-0.8 is drawn in. Similar diagrams can 
he used as aids for the following examples, too. The portions of 
the Mach wave going out from K really require a special computation 
since the average values of X or p for them do not agree, in 
general, with the values of table I. The small deviation was 
tolerable, however. 


9 . FLOWS WITH CONSTMT ENTROPY 3H A PIPE OF VARIABLE CROSS SECTION 


If the cross section of the pipe is not constant, the right side 
of equations (l4) from which it is necessary to start out, here too, 
are preserved. With that, there is the possibility of undertaking 
that integration along the Mach waves which led to equations (20) . 
Nevertheless, the introduction of X and p still remains useful. 
Setting 

a_ ( J_ + JL Ql£\ = M (31) 

a o \ a o w a o & /• 

then 


SA. = -a M (32a) 

fit 

is obtained as the consistency condi ticai for Mach wav9 1 and 

=-a 0 M (32b) 

for Mach vave 2. 

The consistency conditions in the form of (32) contain at any 
given time the differential of only one of the unknown quantities 
X or p while the differentials of both p and v appear in (l^) 
already. This implies an appreciable improvement in the numerical 
calculation . 

The construction of the flow rests op the fact that equations (32) 
are considered different equations. Let G^ be the value which a 

quantity G assumes at the point A, difference G^.- G^ 

and Gjj^ an average value of G taken between A and B. 
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Applying equations (32) to compute from two known points A and B 
the phase at a third, C, which is on the same Mach wave, in the 
form 


X c - \ A « AXc,A = “ M mAC a o At GA 
h c - t^B = m c,B = _M mBC a o At CB 




(33) 


For the determination of the flow Xq and i_iq have to be computed 
by means of these last equations and, at the same time, the position 
ascertained of the points sought in the yt -diagram by the use of 
equations (10) . The calculation process might be explained by an 
example . 


The flow is considered as given along a curve of the yt-diagram 
and, admittedly by X and p (fig* 5, table II). In addition, the 
pine cross section must be s known function of y and t. For that 
it is only necessary to require that F can be differentiated 
with respect to -position and time, a premise which is always fulfilled 
in practice. For this example F is taken in the form 


F * F 0 y^t 


From (31) for M 

M . *.( JL i. + 

a o \ a o y a o7 

The positions y = 0 and t = 0 for which M goes to infinity do 
not belong to this region of flow where such singularities appear 
(for example at the center of spherical waves); it is necessary to 
make special investigations which cannot be entered into in the 
present report 7 . 

The best way to follow the calculation is by means of the 
systematic calculation in table II, To facilitate comparison with 
the description the columns are numbered. The first column contains 
the designation of the point which is to be computed, the second 
column gives the known point which, in caramon with the point to be 
ocmputed, is on Mach wave 1* Column 3 contains the corresponding 

7 Compare G. Guderley, "Starke kugelige Oder zylindrische 
VerdichtungsstSsse in der Mhe des Kugelmittelpunktes oder der 
Zylinderachse Luftfahrtf orschung , Bd. 19 (I9^2),pp.302-312. Tltf-s 
concerns itself with a complicated special case of suoh a singularity. 
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point for Mach wave 2 . The f irst five rows reproduce the initial 
values as well as some further values that hold at the given points 
which are necessary for a later calculation* The calculation of a 
new point is carried out in the form of an iteration method; as 
an example the point It- will he explained. Next, the values for 
Xjj. and (Hj. are estimated (colurms 4 and 5) • Ih order not to 
use too favorable an estimate, it is assumed that X^ = X]_ 

and = (ip. The quantity (X + p)^ is determined for these 
magnitudes and from that, with the aid of figure 3 


farther on 


(A 

V «W 


(f\ 

\ a o/l 


and, 


\ 


and 


(columns 6 and 8) • With these values j "Z. a 

°/4 

are computed (columns 9 and 10 ) « 


\ A 

Wow the average 


\ Q o A 

directions for Mach waves 


1 and 2 f UL 


and 


( 


v - a 


and Ap^ 2 
obtained l 


° /ml, 4 N ° /m2,4 
are formed (columns 14 and 19) and the Mach waves are plotted on 
the yt -diagram. From this y^ and a 0 t^- (columns 12 and 13 ) are 
obtained. With these values (column 11 ) and the average 

values M^i ^ and ^ (columns 15 and 20 ) are computed. To 

continue for Mach waves 1 and 2 Aa Q t^ ^ = a Q t^ - a 0 t^ . 

and Aa Q t^ g = a oA “ a o^2 have to be computed (columns 1 6 and 21) 
and can be* substituted in equations ( 33 ) • The quantities AX^ -j_ 

as well as X^ and (columns 17, 18, 22 , 23) are 

If the values \ and p calculated in this manner do 
nofc #gree well enough with the original estimate, the calculation 
must be repeated in which X and p Just calculated appear in 
place of the earlier estimates. Naturally, the Mach waves must be 
plotted over again, too, in the yt-diagram for this . These figures 
only show the final form at any instant. For that reason all the 
steps in the iteration method are put in the tables. A good view 
of tiie results of the calculation as well as insight into the 
estimates to be carried out by the Iteration method is obtained. If 
the flow is followed simultaneously in a Xp -diagram, as well as the 
yt-diagram (fig* 5, right) . There the X-axis was selected 
slanting up to the right at 45° and the p-axis downward at 45° • With 
a suitable vertical scale X - p^ and therefore v/a OJ is obtained 
immediately on a horizontal scale X + p or W and with the use 
of unequal distributions a/ a Q and P and, for isentropic flows 

p/p too. The X- and p-axes were inclined 4-5° to obtain the 
o 


quantities of physical interest v/a Q , a/a Q , 
system with the conventional arrangement. 


etc . in a coordinate 
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10. FLOW OF AN IDEAL GAS' WITH ENTROPY DIFFERENCES 


The Introduction of X and |i with, the object of obtaining 
equations in only one unknown, at any time, with the iteration 
method for the determination of. Idle flow was possible up until now 

because the expression with constant entropy might have been 

pa 

considered as the differential of a function W independently of 
the characteristics of the incident gas. Naturally, that is no 
longer possible with variable entropy. The computation of the 
flow must, in general, therefore, return to (l4) . The ideal gases 
constitute an exception. Here, as recognized in (30), the function 

W which essentially agrees with / 3-2. for constant entropy, 

d pa 

depends on the temperature alone, and no longer an the entropy. 

If the expression is considered, therefore, in the case of 

pa 

variable entropy as dependent an the variables T and b the 
effect of change in entropy is separated, then the rest 'can be 
written here as a' differential and X and- u can be introduced as 
previously. The change in the entropy along the Mach waves must 
naturally be regarded separately . This is possible without especial 
diff icultles since the entropy is constant along the life lines . The 
transformations are carried through in the following manner. From 
the second law 


T ds = di - i dp 

p 

taking into account (26) and (29) 


1 a:o _ di T ds 

pa a a 



V 


A 


T> 


£ ds 
R 


Introducing W, X, and p as before, the consistency conditions 
are obtained In the form 



dlriF ] + Av T 1 _ ds 

^ / V C p E 6 0 dt 
* 

ding ] + /°v T 1 ds 

dt / . * c p E a Q at 


(34a) for Mach wave 1 


(34b) for Mach wave 2 
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The differential quotients ds/dt formed along the Mach waves 
interfere: The following transformations are possible. Analogous 

to the flow f unction of two-dimensional stationary flows, a 
function \{f is introduced, is constant along the life, lines. 
Uiis can he achieved by requiring that 

U— iJLr (35a) 

H = §- -2- (35b) 

p o p 0 

Along any curve of the yt -diagram 


dk e 3^ + ajf at 

dy ciy dy 


(36) 


Along a life line ay _ y therefore 

dt 


# JL F p v 

& = F o p q " F o e D * 


that is is actually, constant along the life line . At each point 
» of the yt-diagram itself can. be defined by a line Integral that 
leads from, a fixed point A at which might be zero to B • 


\jr 

B 



1 / JL dy - JL ^_vdt 
< - /a V F 0 Pq F 0 Po j 


The physical significance of i{f can be recognized as follows: 
Let C (fig. 6) be the intersection point of the life line through 
A with the line t = constant through B. To begin with, the path 
of integration is along the life line from A to C and, from 
there, out along the line t ** constant to B. Along the life 
line AC, is constant 

+C - - 0 
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along the section CB, dt is equal to zero, accordingly 


...B 


* B -/ ff-u 

B Uc F o °o 


From this, it is evident that represents the mass which is 
enclosed "between the particles at an instant in time for which 
is zero. 


The fact that s is constant along a life line can he written 
with the use of in the form 

s = efllO (38) 

For ds/dt then 


ds _ ds a# 
dt SF dt 


for which £2 is to "be taken, just as ds/dt previously, along 
dt 

the Mach wave considered. 

From (35) and (10) 


4*. JL JLa 
« S'o p o 


for Mach wave 1 


— = a for Mach wave 2 

dt 

o 


Substituting these in equations (?,4), allowing for (23), (28), and 

(29) replacing — -. ft according to (21Tb) by — and p n /p_. 

dm tv d\Jr u 0 

C„. 


By 


u~r 

p o 


yields the following consistency conditions : 
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For Mach vave 1 


dX _ a 
it ■ a = 


/ _v_ 82 nF + _ 1 _ 8 Zn][\ 
\a Q ~Sy a Q dt / 


v o Xpcls 
°Po F 0 


(39a) 


For Mach ware 2 


Si = a 
dt ° 


:'-JL /JL 52 ng + _1_ 52nF^ + T o _F p cbt 
'•» a o v 


O \ a o 


5y . 


St”/ 


C Po *0 




(39b) 


Here P is a function of X + p (fig. 3 ), F/F 0 is Isnown to be 
a function of y and t. From ( 38 ) and ( 27 b) it follows that 


and from this 


it = rt(^) 


d£. w 35L 010 •= constant for a life line (40) 

d\j/ di|f 


For the sake of compactness, introducing 

C v 

T o F _ dk 
N = ~ — — — P r-r 
C F d\If 
0 V 

Then ( 39 ) goes over into the form 


dX * a (-M + N) 
dt 0 


(4la) for Mach ware 1 


$r - a o(‘ M - ») 


dt 


(4lb) for Mach wave 2 


Equations (40) and (4l) supplant the -previous consistency conditions 
(15) and (14) . 
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Before starting the characteristic construction, the problem 
arises here, too, of computing X and p, and now djj. besides, 

ay 

from the initial values. Along a curve K of the yt-diagram let 
the velocity be given by v/a Q , the phase of the gas by p/p D 

and T. From' T with the aid of figure 3 X + p is obtained, 
from v/a Q , X - p; with this X and p are known. Since p/p 0 

are given, and P as a function of T is to be gathered from 
figure 2, « is obtained immediately from (28). As a result of 

plotting it against the values of y from the curve K and 

differentiating is obtained. From (36) and (35) together with 
ay 

(23) g— for the curve K may be computed for the curve K and, 
finally, with that 


dit _ drt &{/ 
dy dy 


is determined. In many cases these computations are superfluous j 
if entrouy differences arise from compression Bhocks, the 

determination of 3* X and p includes their calculation. The 
d¥ 

way the computation of flow hae to be carried out is shown in 
figure 7 and table III with points 4, 5, 6 as examples. The 
related Xp -diagram Is right center. (The points included, in 
addition, in the table and the figures relate to a later section.) 

Along the curve K (points 1-3) X, p, and are assumed as 

d\J/ 

known, in the auxiliary diagram Slj. has been reproduced as a 

function of’ y . The computation of a new point - take point 4 as 
an example - begins, here too, with an estimate of X end |i 
(table III, columns 4 and 5). After that, as before, the following 


are 


computed (X + p)^; (a/a Q ) 4 j (v/a 0 > 4 ; ( 2L±_a\ • f ; 

V a o /% \ a o 


v + ai 


v - ai 


1 a 1 f 1 —z 1 j (columns 6-10, 19 and 24), the position 

V o /ml, 4 \ 0 / m2, 4 

of 4 is indicated in the yt-diagram and y^ and a Q tj + in the 

table (columns 11 and 12) assumed. The determination of with 

dj/ 

the aid of the life lines enters in as something new. It should be 
sufficient for this to draw in a multitude of life lines, simultaneous 
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with, .the construction of the Mach waves and going hack over, these 

to learn the desired value — from the auxiliary diagram. The 

d\{/ 

Xu -diagram is useful for a quick determination of the direction of 
the life lines . The position of the intersection points of the 
life lines with the Mach waves may he estimated there without 
difficulty, and then the average velocity learned. (Compare points Inl- 
and 15 in the yt- and in the Xp -diagram.) After £5. has been found 

d\(r 

and, in addition, P has been learned from diagram 3 (columns 13 
and 14), Mj^ and N^ as well as (-M - N)^ and (-M + N)^ may 

he compu + ed (columns 15-18), the average values (-M - H ^ and 
(-M + N)^ i). for the Mach waves he formed (columns 20 and £5) and 
with ■‘•he aid of A a Q t (columns 21 and 26) from equations (4l) 

compute AX and Ap and, ultimately with that X and p. 

(Columns 22, 23, 27, and 28.) Where the original estimates were too 
had, "t-he computation *-as repeated. 


11. THE GENERALIZED FORM OF THE CHARACTERISTICS METHOD 


An outline shall he given of how to proceed if the simplifications 
given above are no longer possible or if the flew is so small that 
the prepared computations as given at the end of section J do not pay. 
As an example, let the computation of the point 4 he carried through 
from the points 1 and 2 of figure 7. (See fig. 8.) The quantities 

Pl/p o =t 1*44; 7^-= 1.2; T l/ a o ■ °*^ 2 5 

P 2 /p 0 “ 1*866; jr 2 = 1.332; v 2 /a 0 “ 0.400 

correspond to the initial values assumed there. For the medium to 
he investigated p and a must he given as functions of p and it . 

In this case P is obtained, first of .all, from (28) and from that 
and figure 2(h), t. Then 
end (29) . Honce 

e i/ a o = i*02i; 


ag/a 0 = 1.037; 


p/p 0 and a/a o are obtained vith (23) 
P3/P0 - 3-3755 


;; flJL*) - 1.446 

\ 0 ,''l 


P 2 /p o = 1.710; ( 2-r-SL J = 0.637 


a o '2 
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besides 


M 1 = 1 . 747 ; Mg = 1.442 

can be computed. Here, too, an estimate is made in computing a new 
point . For example 


+• «g 

P 4 /P 0 = Pi/P 0 53 3- -44; r k / % = v x /a o - 0.425; = — g 1.266 


With this 


P 4 “ (V^ojK 


1.137 is obtained, whence 


Continuing further 

V a o " 1 - MA » % /o o 


T. = 282.5 
4 


1.390J 



* 1.439 




-0.589; 


(T + a) ml,!|./ a o * 


-0 .613 


With that -the position of point 4 in the yt-diagram may be found, 
giving 


y^ » 1 . 446 ; a 0 t lf = 1.258; A(a 0 t)^^ = 0 . 048 ; A(a 0 t)i^ 2 - 0.092 

and, after further calculation 

- 1.403 

The average values are found to be 

- 1 - 383i ( a/a °L,t ■ 1 - ol75) w * iA » 

KL,* - 1 ' 550i KL A ■ X - oe6i •fe*- 1 '* 3 
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Considering equations (l4) as difference equations then 




+ 






? o (°o\ (to) (?k M a f\ Vj 

p o a o V 0 4,8,4 ' a 4 eA Po P °/ V 


“ a °( 


a-t^ 


" a o t l. 


l)\ 2,1 


Eeplacing 



from (29) and (23) gives 



j°/ £° \ / f2. \ ^1 

a o = C p 0 V p 4l> a 4l,4 p o 

“ ~ %1, h( & o% - a o t l) • 


\ = Jo/ fo\ ( V\ P2 

a o C P 0 ^ p Mn2,lp a -4 i 2 ? 4 p o 


V 2 / v 

— - M 0 1 ( a t, - a t 0 ) 

a o 112 0 ^ 0 2 / 


7 For ideal gases the first term of the left side of (14) may 

he vritten i. ^ISLa, then 0/0 does not have to he computed 
fc dt ’ ' o 

separately. To permit the procedure to he applicable in more general 
cases, this simplification is not used here. 
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Putting in numerical values gives as a result 


0*547 P 4 /p o + v^/a Q = 1.1430 
0.484 p 4 /p q - v 4 /a o = 0.3715 
p 4 /p Q - 1*470 


v. /a = 0.342 
4' o 

Prom the velocity computed above v^/a Q and the velocity at 

a point 4’, estimated for the present, of the connecting line 1.2, 
the average direction of the life line passing through 4 is obtained- 
by an approximation method. If this is proceeding from 4 backwards, 
■Ole more accurate position of 4' is obtained. By interpolation 


between 3. and 2 it^ 1 = = 1.243 is obtained. Since the values 

— , — , do not agree sufficiently well yet with the originally 

j?o a o 

estimated values, the computation must be repeated with the magnitudes 
just obtained as starting values . This gives 


p^/p^ = 1.478; v^/a o = 0 . 3383 ; = 1-243 


32. SIMPLIFICATIONS FOE IDEAL GASES WITH CONSTANT SPECIFIC HEATS 


Generally the flowing medium is an ideal gas with constant 
specific heat or at least can be considered as such, as an approxi- 
mation. In such a case appreciable simplifications are possible. 
Let 


k = C /C 
p' v 


then 


C 


P 


k 


k - 1 




°v 


E 


k - 1 
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From equations (27a), (29),. and (30) 




W = 1 + M- 
2 


/J- 
k - lva Q 


l ) 


With this, it follows that 



X - 2 /-®_ - l\ + 

V 

(42a) 


k - l( a 0 y 

0 



Y 

(42b) 


^ k - 1 f a n ) 1 

a 



0 


and from that 

O ^ 


(42c) 

consequently. 





p=il + ^,1(X + n)| 

Jk-1 

(42d) 

The directions 

of the characteristics are obtained from (9) 

and (10) 


in the form. 


^ = a Q for the life lines 

* fof 1 + k ^ for Mach waves 1 

^ = a 0 ^-l- 4 + 3 -- ^ -- xj for Mach waves 2 
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The consistency conditions for the Mach waves remain unchanged in 
the form (40) and (4l) . M and W are expressed as follows, now. 


M = 1 + 


k - 1 


(X + n) 


X - p 5lnF + 1 

^ T a f 


W 


1 L. 

% 



'4) 


2k 

k-1 


a £ 

d\|/ 


The directions of the characteristics may now he found very 
conveniently graphically. A construction which is suitable if the 
simultaneous treatment of the flow in a Xp -diagram is avoided is 
the contribution of Adam Schmidt. (See fig -90 For the determination 
of the direction dy/dt for a life line, two vertical scales at a 

distance of 1 apart are used with ^ plotted cn the right one 

and & 031 the left one as above . A life line for a phase which 


is given by X and p has the direction of the connecting line 
of the points concerned on the function scales. Similarly, there 


are scales to use for a Mach wave 1, which give p on the 


left and 


1 + “ x 


on the right. For Mach wave 2 




O «. V 

has been plotted on the left and -1 + ■ ■ ^ X an the right-. In 

figure 9, the direction of Mach waves 1 and the life line is given 
for X = 1 .1 and p - 0 .6 . 


If the phases in the course of the construction of a Xp -diagram 
are followed up, the following method~is suitable (fig. 10, right)* 
A vertical line is sent through the 0 -point of the Xp-system and 
the poles P^, Pp, and P^ are determined, where P^ is on a 

level with the origin of the Xp-system and /2 away from it. P^ 
and Pg are directly below and above Pl, respoctively, and 

likewise the distance '/*2 from it. To find the direction of the 
characteristics for a- given phase, a horizontal ray and two rays 


slanting upward and downward at an angle arc tan 


k-1 

2 


are drawn. 


These intersect the vertical line through the origin of -the Xp-system 
at the points Qg, and Q^. The connecting lines P^Q^ 

PgQg, and P^Q^ are the directions of Mach waves 1 and 2 and the 
life line. In figure 10 the construction for point 4 is carried out. 
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This c ms tract ion is especially convenient with a triangle having 

an angle arc tan . Figure 10 and table 4 give an example of 

an application for the same initial values as in figure 7 and with 
Op/°v = constant =1.4. 


13. BOUNDARY CONDITIONS 


If the flowing gas column is not infinite, the variation of - 
the flow is determined "by the phase at the start, in addition, 
also by conditions at its boundaries. For example, a gas can be 
closed off by a piston or rigid wall, flow out into a space with 
a given pressure, or be sucked out of the same. Generally, the 
boundary conditions may be formulated so that relations between the 
phase magnitudes of the gas and its velocity along a curve of the 
yt -diagram ' are prescribed. The number of conditions Which are 
needed for the boundary curve corresponds to the number of charac- 
teristics -which run out from there into the interior of the flow. 
For example, the gas flows out of the end of the pipe into a space 
with constant pressure, with v < a, then the line y = constant 
is the curve for the pipe for which the boundary conditions are 
given. A family of Mach waves spreads out from it inward, while the 
other family and the life lines reach this curve, approaching it 
from within. Hi this case the condition can be prescribed that 
the pressure in the exit section be equal to the outside pressure. 

If the gas is sucked in from outside, Mach waves of the one family 
proceed from the curve of the boundary conditions as well as the 
life lines. Accordingly, two conditions must be given. The one 
states that the entropy of the entering particle is the same as 
the entropy in the outer space, as a second it would be required 
perhaps that the phase of the gas in the entrance section be related 
to the phase in the outer space through Bernoulli's equation 8 . 

(An exact formulation is difficult, since the flow at this location 
is no longer one -dimensional.) If the characteristics of ai 1 three 
families of a given curve lead out into the interior of the region 
to be computed, there are three conditions to prescribe; this is 
the initial value problem already treated. The other extreme, that 
at the boundary of the region of interest, generally, no condition 
can be fulfilled, is physically conceivable, too. For example. If a 
gas with v> a flows in a space at constant pressure, generally no 
characteristic goes inward from the outflow section. 


Compare Schultz -Grunow, loc. cit. 



NACA TM No. 1196 


34 


Actually here — disregarding "boundary .conditions which force 
compression shoots — no effeot on the flow variation in the interior 
is possible from outside. 

The treatment of boundary conditions is explained with two 
examples which are connected with the flow in figure 7. The com- 
putation is entered in table III, aB far as possible. The first 
example includes points 7 to 9 and, admittedly, it has been ass ume d 
that the gas column is bounded by a piston whose life line is 
represented in the yt-diagram as the curve 3, 7, 9. (Whether it is 
practicable to realize such a piston in a tube of variable cross 
section is unimportant for carrying out the computation.) The 
X|a— diagram referred to is in figure 7, upper right. To begin 
with, an estimate is made of the phase at 7 which has been chosen 
^7 = ^3 = 0.800, |iy = U3 ® 0.050. Since the line 3.7 is the life 

of a particle, is already known and is equal to ( — \ , 

W/7 W/ 3 

With this idle values in columns 6-10 and 19 are calculated. As a 
result of drawing in the Mach wave 5 . 7 , yy and a 0 ty (columns 11 

and 12) are obtained and besides vy/a Q from the direction of the 

life line at point 7 which has been reaohed, (This quantity is 
found in column 8 under the value computed from the initial estimates.) 
•Now the quantities in columns 14 to 18 and 20 to 23 may be .computed, 
the value Vy/a 0 obtained from the boundary conditions will be used. 

With that X^ is already known. The quantity (iy is obtained from 
the relation 



Inserting numbers 


0.323 - 1/2. ' (0.417 ~ Hy)j Uy = -0.229. 

Since the first estimate was too poor, the computation must be 
repeated. 

Point S is computed from 6 and 7 by the method explained 
in seotion 9. Prom 8, point 9 is obtained in the way Just described. 

This method of calculation is useful for any laws of motion 
of the pipe; a special argument is necessary only if a discontinuity 
appears. The discontinuity in the velooity is to be considered 
attained on transition of the boundary from a continuous velocity 
variation at very large acceleration. In the yt— diagram that moona 
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that the life line of the piston ■which has a "bend at the instant 
of the Telocity discontinuity is rounded off immediately . Then the 
flow may he drawn accurately just as previously . To obtain 
sufficient accuracy, enough points must he taken on the rounding off 
so tha.t the velocity of the piston does not change excessively from 
point to point, and at each point a Mach wave of the. first family 
may converge and a Mach wave of the second family may diverge from 
there. First of all X must he computed for the converging Mach 
wave and then from X and the Velocity at the incident point ' |i 
determined for each Mach wave. If the rounding off becomes smaller 
and smaller, these points on the rounding off draw closer and 
closer. With that the values of X approach a single value, which 
may he computed from the field before the bend. The Mach waves 2 
spread out in the shape of a fan from the bend and the fan includes 
all values of ,u which lie between the values of p. for the 
velocity before and after the velocity discontinuity. 


For the second example, there is at the position y * y^ an 

open pipe end, through which gas is sucked in from outside and' for 
which two conditions must bo specified along the boundary -condition 
curve. The curve is the curve 1, 10, 13 in figure 7 • In. the outer 
space let 1 t * it-i j for the entering particle therefore drt_ = 0 . 

dty ' 

This Is one boundary condition. As the second boundary .condition ' 
thore is the requirement that the phase in the inflow section be 
related to the phase In the outer space by the Bernoulli equation. 
This condition may be satisfied, already, at point 1, accordingly 

i + v^/2 = i-j_ + v^/2 


or also 



To determine these constants, from figure 3 the temperature T^ Is 
taken foi* (X + [i)^ from figure 2(a) for’ T^, (ii - t 0 )fe^. 




JL_ 

2 


V 

a, 


•.2 


0.292 


Then 
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Sine© (i - i 0 )/a Q ^ is a function of ' T and, therefor©, of \ + n, 

JL = . this boundary condition caii be plotted as the curve K 

a o ^ , . , ' •• 

in the Xii-diagram'- (fig. 7,- lower right). At .bes.t, the computation 
of point 10 begin© anew with an .estimate for and u so that 
the boundary conditions are already satisfied . (columns k and 5) • 
With this, the quantities in columns 6, "7, 8,' 10, and’ are computed 
and Mach wave 2 drawn' in with’ that . The quantity a^t' 10 iB 

obtained in column 12. the -values y lf) - y-i . and - M ss o. (Columns 12 

.. ■ ■ ' 

and 13 are given beforehand . ) Now the quantities 'in cblumns 14 to 18 
can be obtained. 

To determine, with this, -the quantity (-M +. N.) in column 25 
it is to be noted that (-M + N) • for the particle originally in 
the pipe has the 'value, perhaps,, at point 4 and changes die- 
continuously for the particle recently sucked' into the quantity 
(-M + N) 10 . 

On that account the life line is drawjl, which separates the 
particles in the interior originally from those particles flowing 
in from outside. This intersects!” Mach wave 4,.. 10 at point 11. 

Then the following- is ’obtained’ (column : 25)‘ 

t 

tA(a 0 t) llil0 (-M + N)^ 

The quantities in columns 26, 27, and 28 may be computed now. As a 
result of inspecting the curve of the boundary condition in the 
X|i“diagram with the value of u ' fotnd, "K is obtained (column 23) . 
“Hie computation is repeated with the values found in this way. 

From points 6 and 10, point 12 is obtained in the manner 
described in section 9. In connection with that the difficulty just 
described appears again in finding the average value for (-M + N) . 
From 12 and the boundary condition’, point 13 may be computed by 
the method just presented. 

The \(i- diagrams of the two last examples were kept separate 
from the >41 -diagram drawn for points 1-6 for the sake of clarity . 
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If the various figures are visualized as being joined - the upper 
diagram' connected to 'the middle one at the line 6, 5, 3, the middle 
one with the lover one at the line 1, 4, 6 — it is recognized that 
the plane is covered with several sheets which are connected ai opg 
the figures of the characteristics . There is such a superposition, 
already, in the lower Ip. -die gram; there are to he imagined inclosed 
the quadrilateral 10, 4, 6, 12 along 10, 4, the triangle 1, 4, 10, 
along 10, 12 the triangle 10, 3, 12. 

In addition to the boundary conditions, transitional conditions 
can also appear in the interior of the flow. In the example just 
discussed just that would have been the case, if in the outer 
space it were different from rt, . At the location of such a 

discontinuity for it agreement of pressure and velocity must be 
required. To go into such questions with greater detail lies 
beyond the scope of this report. 


14. TRANSITIONAL CCNDITICNS AT COMPRESSION SHOCKS 


The flow in a given part of the yt-plane is defined by the 
initial and boundary conditions and is calculable by the methods 
derived up until now. It is possible that it might happen during 
the coinstruction that regions of the yt-plane ore covered with 
phase quantities several times . This is the sign for the appear- 
ance of compression shocks . The entropy is no longer constant 
after the passage of a compression shock. On that account the 
computation of compression shocks simultaneously includes the 


determination of the function s(\j/) or ~(v|/), too, for the 

dy 

region of the yt-plane behind the compression shock. 


For the mathematical treatment, a compression shock is to be 
considered a curve along which two flows collide, which are related 
to one another and to the diroction of this curve by transition' 
conditions. It will be the ^problem of this section to derive 
these (known of themselves) " transition conditions in a convenient 
form, for the present purpose . 

Proceeding from a stationary compression shock, that is from 
a compression shock whose front is at rest relative to the coordinate 
system selected, let the index I designate the phase before the 

9 Compare Ackeret for instcnco. Beitrag Gasdynamik in Handbuch 
der, Physik, Bd. VII, p. 324 and following pages, Berlin I927. 
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shock, the index II the phase after the shock. The additional 
index s might point out that this concerns the calculation of 
a stationary shock. Then the momentum, and the energy theorems 
as well as the equation of continuity are written in the form 


hs * °I8 T IB 2 " *118 + toW 


(43a) 



+ I v ^ 
2 3# 


IIB 


+ 1 v 2 

2 IBs 


(43b) 


O v = P V 
Is Is Ils IlQ 


(43c) 


Furthermore, the characteristics of the gas concerned must be known, 
possibly in the form 


P - p(i, o) 


(43d) 


If the quantities in advance of the shock ij s , Oj s , and Vj 0 are 

known, then the compression shock is therewith calculable , Actually 

all three quantities enter into the general gas. laws, too,- as 

parameters. In order to carry out the computation practically, in 

such a case, p T _ from (43c) and i_ from (43b) have to be expressed 

1X3 IIS 

as functions of vjj s and the known quantities and then substituted 

in (43&) • With that, an account of (43d), t> , too, is a 

‘ Ils 


function of v. 


Ils 


and the known quantities in advance of the shock. 


In this manner an equation for v. 


Ils 


alone is obtained which must 


be solved numerically in a suitable manner. For an ideal gas for 
which c is not constant, equations (43) transform with the aid 

of (23, as follows : 


Is 

3 IIs 


ET. 


Is 


+ £i§_ Y 


Ils 


Is' 


2 = 


!.( T I8> + tlB 


v Is 



ET II8 + ’ns 2 

(44a) 

i ( T IXs) + J T ile 2 

(44b) 

T IIs 

(44c) 


P 
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Since o appears here only in the combination Pjjg/Pjg ^^-7 

Tj g and v^. g still remain as parameters upon which the phases 

behind the shock depend. To calculate the shock curves numerically, 
it is useful, first to regard Tj s and T-j-j b as parameters and 

determine Vj a from this subsequently . The computation process is 

the following: Prom (44a) and (41(c) . - 


RT 


I 

Is 


HTit 

+ V T = =44 + V TT 

V x „ I 0 Vjj I ^- 3 


(45a) 


As a result of squaring this 


R 2 T- 2 


2 2 

I f- + 2 RT Ib + v Is 2 = + 2RTxx„ + v-, 2 


’Is 


'lie 


L IIs v IIs 


(45b) 


Introducing 


" i IIs ‘ 1 Is 


gives 



V - ^ 


(W) 


from (44b) . 

Putting this in (45b), the desired equation for -IS 2 is 
obtained as 


Is 


^ j^2Ai - 2R^T' tt _ - T To ) I + v To 2 I -4Ai 2 + 4R£ifc 


ns is, 


a 


Is 


=) 


■ s2 ( T IIS 2 - T ls 2 


e ) 


l ns ■‘•is 


-2R 2 T Is 2 Ai = 0 


If Vj s is determined, then v^^. g and Pjjg/Pjs ,are cani P u ‘* :i ® < i ^ 
turn, with the aid of (46) and (44c) ; finally 
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For an ideal gas yith constant specific heats, the following trans- 
formations may be undertaken . According to the familiar relations 


and 



a? - kRT 

Equations (45a) and (44b) are written in the form 


^s 2 + 
^Is 


a TTn 2 

T Is = £ 7 - + T IIe 
^IIs 


k - 


T G ls 2 + T le 2 


~ a IIS + V IIS C 


or 

1 1 , Y Is _ l/ a IIs\ 2 1 ■ + Ills 

k t Ib^®Ib a Is k \ a Is ) T lls' a Is a Is 


(47a) 




(47b) 


By this, aq; Is /a l8 and and, with that, the other quantities, 

too, depend on the parameter v ls /a-j- s alone. 

\2 

To compute v IIs/ a l8> ( a IIs/ a Is) ie eliminated: 




Is obtained as a result. 
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The solution of 'this equation is found, immediately, if it is 
home in mind that on account of the form of (47) a solution is 
represented hy 


r IIs/ a Is = T ls/ a Is 

then 

7 IIs „ a f is. + k - 1 T Is\ 
e Is t+1 \ T lB 2 «!. ) 

Using this, the following is obtained from ( 47 b) 



°ne /p l3 = ( T XB /a is)fls/ T ii s ) 


P IIS /,p lB 



= v_ /a a^ /v 
Is' Is Is' Us 



The change of entropy is of interest, as well; with the aid of (27) 
and (28) , these expressions result 


s - s 

-J I .3 :ia 

E 


id! ina\ - 


*- 1 w 


te£tla 

p Is 


fe-L S Ie „ _gk... 2n!lle - j n IlS. + ln Ills . 2 zn Ills 


B 


k - 1 a -r a-r ■ 

Is Is 


*Is 


a Is 


__ 2 — Zn _ r . II , B « Zn Il£ + zn ZEES. 
k “ 1 a ls a ls a ls 
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From thiB 


rt IlB 
n Is 

Arbitrary compress Ion shocks result from the stationary compression 
shocks Just calculated because a Telocity is superimposed. In doing 
so, the thermodynamic phase quantities before and after the shock 
for •which accordingly the index s can be omitted are retained and 
moreover the velocity differences. Since the phase in advance of 
the shock is already given in the construction of flows, before the 
shock is computed, the relative velocities with respect to the phase 
in advance of the shock are formed. Let 


-( a ils\ " Ils fjEs 
"Ws ) a iB V IIS 


u absolute velooity of shook front 

Au relative velooity of shook front with respect to partioleB 
in advance of shook 


Then 



Av, 


II, I 


= v. 


n 


- V, 


2 / Y Is 

k + 1 1 



The signs appearing in this are not astonishing. A stationary 
compression shock in a gas which moves in the positive direction 
propagates itself in a negative direction relative to the material 
ahead of the shock, and in so doing, produces a change in velocity 
in the direction of its propagation velocity, that is, in the 
negative direction, too. Naturally, compression shocks, which travel 
in the positive direction in the material at rest are also possible, 
the Bigns of the velocities have to be changed for these . The 
thermodynamic phase quantities of this are not touched upon. Corre- 
sponding to the distinction which had been met in Mach waves, these 
last compression shocks are designated compression shocks of the 
first type, those which propagate in the negative direction as 
compression shocks of the second type. In figure 11 the pressure 
ratio, for an ideal gas with k = 1.405 the propagation velocity of 
the compression shock and the change in entropy (expressed by 
«Ij/jti) has been presented as a function of the velocity change 

^ T II I* ^ or com P ression shocks of tiie first type Au and Av^ ^ 

are to be taken with positive sign, for compression shocks of the 
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second type with negative sign. The fundamental numerical values 
appear in table Y . Such a diagram would have to be used to apply 
the characteristics' method in the form given In section 11 in the 
computation of compression shocks. How are these relations for the 
compression shock expressed in terms of X and |i ? If two 
compression shocks which only arise separately from superposition 
of a velocity - they are distinguished by the indexes a and (B - 
are represented in a X,fi- diagram, that is, if the phases In advance 
of the shock Xj a ; n-j. (J-jp and the phases behind the 

shock are plotted, then here, too, the eapression must be arrived 
at that the thermodynamic phases in advance of and behind the 
shock, as well as the velocity differences for both compression 
shocks are the same. Accordingly, 

\a + ^I,a = X IP + 4 I|3 



By subtraction of the first two equations 

\XC,a " X I,a) + (^H,a " |1 I,a) = ( X II,P ’ 
Rearranging terrs in the third equation gives 

^II,a " X I,a) (^II,a “ ^I,a) = (*113 ’ \,&) 

From the last two equations it follows that 

X " X * X * X 
II, a I, a II,p I,p 

X II,a " M I,a " ^11,3 ’ ^1,3 


+ (“ns • “is) 
' ("ns ' "is) 
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that is, the changes In X and }i In a compression shock are main** 
tained In the superposition of a velocity* Accordingly } the shocks 
are designated by 


^11,1 ■ x n " x i 


and 

^11,1 = ^11 " 

The following relations hold for ideal gases with constant specific 
heats, according to (42) 



is to be computed from Xj and m-j by (42c) . For the 
expressions in curved brackets 
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ere Introduced. These quantities, as well as Au/a-j-, itji/rtj, 
and Pjj/pj defend only on v j s / a j according to relations 

previously developed. They are plotted in figures 12(a) and 12(h), 
and, admittedly, the upper designations refer to tlie compression 
shocks of the first type, and the lower designations to compression 
shocks of the second type. Figure 12(h) represents an increased 
section of figure 12(a), with ,the. apnronriate numerical, values in 
table V. ' ' 

The following example shows a first application of this diagram. 
Xn a pipe of constant cross section there is a quiescent gas of 
constant entropy and constant pressure, the sonic velocity is taken 
to he a^j. = a Q . Suddenly, a piston is driven into the pipe at a 

uniform speed of 0 . 5a Q . What is the ensuing flow like? Figure 13 

shows the yt-diagram. • The starting point of the piston motion lies 
at the origin of the coordinate system. The life line of the piston 
is shown with hatching. A compression shock forms in front of the 
piston, which imparts th.6 velocity of the piston to the particles, so 
that the particles behind the compression shock move with constant. • 
velocity . Corresponding to the phase in front of the compression 
shock is 


Xj = 0 j Pj — o 

The velocity behind the compression shock is 


therefore. 


T n = 0 -5 a o 


e^n ' ^n) ■ 0,5 


x n ' 1 

From this, on account of X^. = 0 and p^ = 0 


AX - Ap =1 
II, I II, I 


Since a j/ a Q = 1 this gives 
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As a result of causing this straight line in the AXA^-diagram 
(fig. 12(h)) to intersect the shock curve, the following is obtained: 

AX = 1 .022 j Ap - 0 .022 ) & = 1.346) itjj/icj » 0.970 

X I]; = 1.022) p * 0.022) u » 1.346 

From Xjj and Pjj, P is computed by (42d) , from this by (28) 

p n /Pl = 1-970 


The goal would be reached somewhat quicker in this by application of 
diagram 11. 


15 • PRELIMINARY ARGUMENTS IN THE DETERMINATION OF A COMPRESSION 

SHOCK IN THE FLOW FIELD 


It is the object of this section to show first of all by what 
data a compression shock In a flow is determined, and, secondly, to 
give a method by which the computation of such a compression shock is 
possible . 

. Ap can be readily shown, the velocity of a compression shock is 
larger than -the velocity of a Mach wave in the material. This means, 
that the flow field in advance of the co uprose ion shock remains 
unaffected by this and can be computed independent of it. It will be 
assumed to be knovn vhat follows . For the field behind the shock, 
a c empress ion shock of the first type represents on the one hand the 
start of life lines and Mach waves 2, on the other hand the terminal 
of Mach waves 1. It follows, from this, that the flow behind the 
Bhock and the shock itself are mutually related and can only be 
computed together. This Is the reason, therefore, that the computation 
of the compression shocks becomes, essentially, more complicated 
than the computation of other parts of the flow. 

Next will be shown how examples can be conceived of flow fields 
with compression, shocks . If in the yt-diagram (figs. 14(a) and l4(b) ; 
the flow field in front of the compression shocks and the portion CD 
of the life line of the compression shock is given, then the phases 
behind the shock are also determined. From the slope of the life 
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line the propagation Telocity of the compression shock is given, 
namely for each point of CD. Beside, the phases in front of the 
shock can be learned for the points of CD; with this the phases 
behind tho shock are calculable . From the phases behind the shock, 
a portion of the flow field behind the shock, namely the region CED 
(fig. 14(a)) may be computed, or if the entropy is known for the life 
lines at the lower end of C. The region CFD (fig. 14(b)) as well. 
It is necessary to go forward along the life lines and Mach waves 
backwards along Mach waves 1. Imagine in figure 14(a) that the 
computed life line CE is realized through the motion of a piston, 

then there is a flow in which a compression shock appears and which 

satisfies a boundary condition (if not prescribed, too) . In 
figure 14(b) it is necessary to imagine another flow field adjoined 
continuously at the lover end of CF ; here the compression shock 
and the flow determined by it satisfy the condition that it is 
compatible along the Mach wave CF with another flow. 

From these flow fields the following is recognized; the 
compression shock through the portion Cl of the life line of the 
piston or CF of the Mach wave is defined as far as It is reached 
by Mach raves of Its type (here the first, therefore) . A change of 

the life line of the piston outside of CE or the Mach raves 

outside of CF propagates along Mach wave 1 in the yt-diagram, to 
be exact, and neglecting cases in which a second compression shock 
arises, attains the compression shock at the upper end of D, 
certainly. On the other hand a change brought about between C 
and E or between C and F in the boundary or junction 
conditions takes effect at that position an the compression shock 
where the Mach wave 1 concerned reaches it, that is, the portion CD 
is certainly changed. 

If the life line of the piston is known beyond E to G or 
the Mach wave beyond F to H, then a further portion of the flow 
field is thereby determined, without the necessity for knowing the 
continuation of the compression shock beyond D; it concerns the 
regions C1GJD or CFHKD . 

It will now be shown how to procede fundamentally to compute 
a compression shock for specified boundary or junction conditions. 

As a concrete example assume the compression shock to be produced 
by a piston which experiences a sudden jump in velocity. (See fig. 15*) 
The starting point of the compression Bhock is that point of the 
life line of the piston at which the velocity jump appears. The 
phase immediately behind M can be ascertained immediately by the 
method applied to the example of the last section. The compression 
shock - as in previous examples of Mach waves - is computed in 
individual sections, which are so small that the phase quantities 
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for them, may "be regarded as varying linearly. Ab just carried out,, 
the phases behind the compression shock are calculable, if the 
velocity of the shock is known. The velocity at M is known. 

Along the portion of the compression shock to bo computed, M, N, 
the phase change and, with it, the change in propagation velocity 
of the compression shock, too, are regarded as linear. Accordingly, 
for all possible shocks which satisfy the transition conditions, 
the portion M, N, of the compression shock depends only on a 
single parameter, the velocity change between M and N, to be 
exact. As a result of compu+ing the field behind the compression 
shock for various values of this parameter, by interpolation, that 
shock may be ascertained which is consistent with the specified 
piston movement. At best, for this N is permitted to travel on a 
fixed life line in the field in advance of the shock. Let C be 
the point on the life line for which the Mach wave 1 passing 
through N proceeds. Now the region OPQN may be computed in a 
familiar manner. For the determination of the extension of the 
compression shock NR the phase behind the compression Bhock at 
the point N may. be regarded as given everywhere along the entire 
Mach wave NQ. . On the other hand, that value of velocity changes 
between N and R . has to be determined by interpolation, which 
relates to a flow field that continuously joins the known field 
along NQ. 

With these two typos, namely the computation of a compression 
shock going out from a piston or wall and the computation of a 
compression shock continuing into or arising in the interior of the 
flow, the most impor tan t problems have been mastered that can appear 
here. The Interpolation methods described become pretty tedious; 
instead of them, iteration methods will be used, which actually lead 
to tho goal more quickly. The interpolation method was mentioned 
previously, however, since it affords better insight into the basic 
relations . 


l6. EXAMPLES OF THE COMPUTATION OF COMPRESSION SHOCKS IN THE 

FLOW FIELD 


Example s trill be given of how the problems formulated in the 
preceding section can be solved by means of iteration methods. Let 
the flow be that, computed in figure 10 and table IV. As the start 
of the new portion of the compression shock to be computed, point 1 
is chosen in every case, accordingly it is identified with the 
point M (fig. 15) once and with the point N a second time. The 
new portion of the compression shock to be computed that corresponds, 
to MN or NR, accordingly, is assumed fco end 6n the life line 8, 9 
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of figure 10. The phases in front of the shock for N or E are 
obtained as a result of interpolation along this line. For these 
calculations it is necessary, on that account, to have the know- 
ledge of the flow field in front of the shock at the points 1 
(M or IT ) and 8 and 9 • In table VI which has the same arrangement 
as table IV these values have been recorded. While it sufficed to 
know for the construction of the flow field, here jt Itself 

at 

must be known. These quantities for points 1, 8, and 9 are located 
in column 26. In the designations, in these examples, the only 
deviation from figure 15 is that only points on the compression 
shock are characterized by letters . Numbers are used for points of 
the flow field, corresponding to previous use. 


We begin with the more elementary problem of continuing a 
compression shock in the interior of the flow. For this the phase 
behind the shock at the point N and the phases along the Mach 
wave Njj, 10 (fig. 16(a)) may be considered known. The phases 


at and at point 10 appear in table VI, phases in between are 

found by linear interpolation; moreover, for N^ the 'velocity of 
the compression shock and it have been given, (columns 25 and-26) • 


Besides 


djt 

d\J/ 


for the life lines lying below N may be viewed as 


computed. It was entered for point 10 in the corresponding column. 
If the distances between points on the compression shock are not 
chosen too large, it iB sufficient to regard between them as 


as constant. In the following this has happened throughout. Since 
N^j and 10 lie on a Mach wave, the consistency condition must 


naturally be satisfied. 


In connection with the flow calculation the existing data are 
to be taken from the preceding calculation steps. The real 
computation begins with the fact that the difference in from 
its value at the starting point of the portion of the compression 
shock to be computed (N here) is ascertained for the life line 
up to which the compression shock is to be computed (8, 9 here) . 
This computation is carried through along the curve of the initial 
values in. figure 10, the life line 8, 9 used here passes through 
point T there. By (37) 
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By (23) and (42d) 



1 + 


+ n) I 


2 

k-l 


Just as for figure 10, F lias the form 

F = F G y 2 t 

For point 7 


y = 1.450; a.t = 1.180; X = 0.66; n » -0.16; rf = it = 0.849 
0 7 o 

For N the corresponding values appear in table VI. With this the 
following is obtained: 




p ° H 


1.110 






0.930 




= 2.425; 



0 




1.020 


Mr = 2.425 x 0.075 + 1.020 x 0.03 = 0.2122 
7,N 
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la figure 10 had already "been given, it must be the same as 

that found from, the quantities just computed. In fact 



= 0»Qkg = 0.230 

Mrj N 0.212 


Shis is the average value of M. as can he gathered for the 

stretch 1.7 from the auxiliary diagram in figure 10. After these 
preparations, the actual iteration method is reached. To begin 
with, the phases at the points and 11 are estimated, in 

that 11 is the intersection point of the Mach wave 1 leading 
backwards from R with the given Mach wave IT , 10 . Since no better 
reference point exists for the estimate, these phases are equated 
to the phase at H_. Moreover, still another estimate is needed 

for — behind the shock; for this, the same va3.ue that prevails 
dijr ' 

at the lower end of F is chosen. With these assumptions, the 
figure F, R, 11 may be drawn in figure l6(a) . Starting with the 
life line of the compression shock KR, whose direction here is 
the same as the direction of the compression shock at F (table "VI), 
R Is obtained as the intersection point with the life line 8, 9 • 
Then the Mach wave R,ll is drawn in proceeding frcau E backwards. 
The direction of this Mach wave was taken in the familiar manner 
from a Ap- diagram (not given here) . From this f igure the position 
of E in advance of the shock is learned by interpolation along 
N,10 the phase a+ 11. (See table VI.) From this may be obtained 
the values entered further on in the respective lines which are 
necessary for later computation. Proceeding from Xjj by means 

of the consistency conditions, the quantity A-^jj is computed for 

the Mach wave (ll/R.^.) . For this the initial estimates for idle 

phase in E TI are taken as a basis and then columns 6 to 13, 17, 

15, 1 6 and 18 to 20 computed. For X ^ so obtained the 

properties of the compression shock are taken from the shock 
diagram 12(b). The following computations are essential to this 


* A„„ 
Eli, I Eli 



0.95^ 




= o.95fc/L.o4o - 0 . 91 & 
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From the shock diagram. 


^ = 0.0130; "eh/xh! = 0 . 9 T &-, 


= 1.307 

®BI 


From this it is computed that 

Au 

Au - 0o013?; = 0.830; = 1.36b 

EXI.T EH a 

’ o 

fSn = '°- 090J U R /a o = 7 SI /a o + A V a 0 * 1 - 661 

A portion of these results are given in table VI (columns 24 to 26). 
Moreover 


dJt „ *BII " "Nil _ 0.830 - 0.781 = 0 230 
47 A\!/ 0 .2122 

T?,N 

To inprove these values, let a second iteration step be carried out. 
First, the figure N,B,11 has to be drawn again for the values 
just obtained. The average direction of the compression shock is 


“tiiH.K = 2 ( % + ^ = 1 ' 733 


Then B^ and 11 are obtained by interpolation , X^jj from the 
consistency condition for the Mach -wave 11,B TT . 

To find the characteristics of the shock, it is necessary to 
carry out the following computation 


AX, 


Eli, I 


= 1.456 - 0.493 = 0.963; AX^ « 0.927 


From the shock diagram 


Au 


A)I = 0.0130; it /*__ = 0.980; —fi - 1.310 
EXX J KX 8^ 


From this is obtained 

(i = -0=087; it = 0.828; 

R,II Eli 


u /a = 1.657; & = 0.220 

El Jj 0 
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An additional iteration step is not necessary any more . In the 
second example (fig. 16(b)) the compression Bhock is produced by 
the sudden Telocity change of a pis ten. The point of the yt- diagram 
at which this velocity Jump takes place - let it be designated M 
in agreement with figure 15 - is to coincide with point 1 of 
figure 10. From the point M the piston has the velocity corre- 
sponding to the life line in the field in front of the shock, in 
particular the velocity at M in front of the velocity Jump is 
0 .425 a 0 • At M the velocity changes, suddenly, to the value 



0«925a o and rises until the instant a Q t 


1.3 to the 


uagnitude O.975a 0 . This and the flow field as determined by the 

initial conditions and. the piston motion up to the point M is given. 
Next, the phase behind the shock at the point M is computed. 


^KEI,I 

a 

o 


\ a o/M,I V a o/ M>II 


0.500 


X 

i(M - An) 

2 

As a result of this line in the shock diagram 12(b) intersecting 
the Bhock curve, the following is obtained 

AX = 0.986; Af I = 0.020 

W> MI = 0 -9T7; £21- 1-333 

“MI 


0.483 
= 0.483 


Vn - 1 ' 620 ’' “m.h - -°' SS9 

”m,h ’ 0-78l < ^-1-805 


From this 
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The phase at is known with that, (table VI.) Now the 

difference must be commuted, over again, from the life line of 
the piston for the life-line up to which it is desired to compute 
the compression shock. It is desired to allow the compression 
shock to end at the life line 8, 9, here too and take the phases 
in 8 and 9 (table VI) from the' preceding example' and 


r - Ap- = 0 .2122 

8,M ^n,m 

The computation of the compression shock makes use of figure M,S, 

N, 11 . (See fig. 16(b)). M, S, N iB the life line of the compression 
shock; N, 11 is the Mach wave 1 returning from N; 11 * S is the 
Mach wave 2 returning from 11 . To begin, an estimate of the phase 


at the points N^, 11 and S TT is made and this is chosen equal 

everywhere to the phase at Mj-j-. In addition, an estimate for 

is necessary. Let ~ -a 0.230 as a start. Figure M, N, 11 , S 

dy 

may he drawn with these assumed values . The order in which the 
points were named corresponds to the order in which they came up 
in Idle drawing. For the positions of N and 11 obtained thereby 
the phase in front of the shock (see table VI) or the velocity of 
the life line is obtained by interpolation. The iteration method 
begins at point 11 and it can be shown that p, can be only 

slightly different from pg because the line element Sjj, 11 is 
Bmall relative to the other’*’ ciimens ions . The quantity p c can 


differ from p. 


M- 


II 


3 II 

only slightly, since it originates in linear 


interpolation between 


Therefore 


^11 = %!! 


M and N', and N lies very close to M. 
is chosen as a starting point. If the 


velocity of the piston at .11 that is known from the boundary 
conditions is used for this may be computed. 

consistency condition for idle Mach wave 11 ,N 

Now the following computation 




From the 
is obtained. 


^11,1 = 1,017 ^N = 0,988 


and from the shock diagram . 

^=0.020) = 0.973; jjk- 1 . 33 * 



NACA m Wo. 1196 


55 


from this 


"11,11 ■ '°- 077! 
7 > ' 


^ I3 . = 0.827; - 1.678 


Further it ie caD.cu3.ated that 


& - 0.217 

The phase at S is obtained by interpolation between M and N . 
With the aid of the consistency condition for the Mach wave S,ll,H 2 i 

is finally obtained, and X^i from the boundary condition for 


point 11. The' first iteration step ends with that. It is necessary 
to check whether the quantities X-, ^II* %H: 

and 




computed agree sufficiently with the original estimates . 


To increase the accuracy a second iteration Btep might be 
carried out. On the basis of the values just computed, the figure 
is redesigned and the computation is carried out in the manner Just 
described. The value for just computed is taken as a beginning 

The following calculation is obtained for the determination of the 
characteristics of the shock 


AX 

HU, I 


1.003; 



0.967 


From the shock diagram 

^ = o.° l8 i ,1 = 0.975; ^ = 1.327 

' ' ®HI 

from this 


-0.083) 0.828) 1 - 1.678) 


ch|/ 


The computation is continued in the manner given until the phase 
at point 11 is obtained, again. An additional iteration step is 
not necessary. 
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IT. SUMMARY 


The differential equation system for nonstationary, one- 
dimensional flows possesses three families of characteristics; the 
thermodynamic and the flow phase are described by three variables. 

As a result of setting up consistency conditions for the charac- 
teristics passing through the point for which the conditions have 
been set up, three equations are obtained from which the phase 
may be obtained. Jn that a possibility for the computation of the 
flow has been given fundamentally.' The report carries out these 
ideas, in general, and brings the simplifications which are possible 
under special assumptions, as well as detailed examples. Compression 
shocks appear., in this, as transitional conditions in the interior 
of the flow and are, likewise investigated in detail. 


Translated by Dave Feingold 
National Advisory Committee 
for Aeronautics 
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Figure l.~ Curvilinear coordinate system | , tj . 
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Figure 9 
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Figure 11.- Characteristics of compression shocks K = 1.405. 
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Figure 14a. 


Figure 14b. 





